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This paper is a treatment on linear and planar phased arrays of current 
sources, whose amplitudes are uniform and scan-invariant. By recognition 
that the radiation impedance of an array element is an analytic function 
of a complex scan variable, a powerful mathematical tool becomes avail- 
able for the investigation of some important properties of the impedance 
as a function of scan. For example, it is proven that in a finite array the 
impedance seen by such a scan-invariant current source cannot be per- 
fectly matched over a continuous scanning range using lossless, linear, 
passive and time-invariant elements. This result is extended to the infinite- 
array case by treating the latter as a periodic structure, and assuming 
that the Green's function of the unit cell is analytic with respect to the 
scan variable. The theory includes both linear and planar arrays. Among 
other results it is shown that the element impedance in an infinite array 
must be of a specific mathematical form. It is hoped that by recognizing 
the limitations imposed thereby, useful guidelines will be established for 
achieving optimal match of an array into space. 

I. INTRODUCTION 

The class of antennas widely known as phased arrays includes es- 
sentially two types of radiators: stationary and steerable ones. The 
first operates at fixed amplitude and fixed relative phase between the 
array elements. Consequently, the antenna characteristics, such as ra- 
diation pattern, input impedance, and mutual coupling between ele- 
ments, remain unchanged during the entire operational lifetime of the 
antenna. The steerable antenna is characterized by time varying ex- 
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citation. The relative phase between adjacent elements is varied either 
mechanically or electronically to bring about a variation in the orien- 
tation of the beam. In most instances scanned arrays are large, in size 
and may contain several thousand elements. Their illumination has a 
linear phase taper. As a result the antenna characteristics become scan 
dependent. The relationship between scan angle and various param- 
eters of interest such as gain, element impedance, and mutual coupling 
between elements have been the subject of intense investigation in re- 
cent years. 1 - 2 One particular direction has been towards improvement 
of the impedance match over wide scanning ranges. 3 At present the 
merit of a matching technique can be determined only relatively to 
other techniques. To the best of the authors' knowledge an absolute 
mathematical criterion, based on physical realizability requirements, 
has not been formulated. Some investigators 4 - 5 claim that a perfect 
match of an infinite array for all scan angles (at which the active im- 
pedance is not infinite, zero or purely reactive) can be achieved by an 
infinite set of interconnecting network elements. However, the proof 
is based on the assumption that the scan-dependent equivalent load 
impedance at the array-space interface remains unchanged after the 
sources have been interconnected by coupling elements. Although this 
assumption has been successfully applied 3 - 5 to improve the matching 
capability of an infinite array, it is incorrect to use it in a perfect 
matching scheme. 

In this paper a new mathematical approach to phased array anal- 
ysis is presented. The model for the analysis is a phased array of ideal 
current sources (electric or magnetic) of scan-invariant uniform am- 
plitude. This model is further discussed in Section II. The analysis 
itself is based on the general laws of antenna theory and on those 
properties which are common to all phased arrays represented by the 
model. 

The first part of the theory is devoted to finite arrays and is treated 
in Section III. The starting point of the theory is a theorem which 
establishes that the radiation impedance of an element in a finite array 
is an analytic function of the scan angle. Further, it is shown that an 
element in a linear or planar phased array cannot be perfectly matched 
over a continuous scanning range by using lossless, linear, passive and 
time-invariant elements. Then it is demonstrated that the directions 
in space of the beams' maxima are eigenvalues of a Laplacian differ- 
ential operator with periodic boundary conditions which are related to 
the phase taper of the array, and several useful properties of those 
eigenvalues are derived. 
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The second part of the theory appears in Section IV and is devoted 
to infinite arrays, which play an important role in the analysis of large 
phased arrays. The investigation is based on a transformation between 
the scan angle and a complex variable s = a + j&, which can be in- 
terpreted on < a ^ 1, £ = as the trigonometric sine function of the 
angle between the plane of the array and the direction in which a chosen 
grating lobe propagates. It is subsequently shown that the element 
impedance, as a function of s, is restricted to a specific mathematical 
form. Recognition of the limitations imposed thereby may provide new 
insight into the behavior of such arrays. 

II. PRELIMINARY REMARKS 

The model chosen for the following treatment is a linear or planar 
phased array excited by a set of ideal current generators of uniform 
amplitude and linear phase taper. The description ideal implies that 
the sources have no internal impedance and are invariant under any 
loading. This means that except for the relative phasing between con- 
tiguous generators the currents are scan independent. Frequently in 
antenna analysis induced currents are replaced by equivalent sources 
by application of the equivalence principle. Such currents are not part 
of the sources. The induced currents are accounted for automatically 
by fulfillment of the requirement that the tangential component of 
the electric field has to vanish on all conductors. In general, the source- 
current amplitude in each clement of the array may be a function of 
scan. However, this dependence is generally unknown and is often 
neglected in theoretical work. The types of excitations commonly used 
are the "free excitation'' and "forced excitation".* The first assumes a 
generator with a scan-invariant internal impedance which is capable of 
delivering scan-invariant incident power. In the latter a constant termi- 
nal voltage or current is maintained. As pointed out by Oliner and 
Malech free excitation is easier to realize in high-frequency technology 
than forced excitation. The latter, however, is more tractable here. 
The results of this study remain valid for scan-dependent excitation 
as well, provided the current density of the source is a smoothly vary- 
ing function of scan angle and can be analytically continued into a 
complex scan-angle plane. 

Under the assumption that the array is excited by a uniform ampli- 
tude and a linear phase taper, the current density excitation function 



* A. A. Oliner and R. G. Malech, Rof., 1, pp. 209-211. 
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of an M-element linear array (Fig. 1) is given by 



J(x, y, *, id = 1 



J (.r - ma, y, z)e im * , ma ^ x ^ (?n + l)a, 
m = 0, 1, ••• , M - 1, 



(1) 



J(a, 2/, 2, ik , ft,) = 



otherwise 

and that of an M X N element planar array of rectangular symmetry 
(Fig. 2) is given by 

(.r — ma, y — no, z)e 

ma ^ x ^ {m + l)a, 

nb ^ y ^ (n + 1)6, ^ 

m = 0, 1, 2, ••• , M - 1, 

n = 0, 1,2, ••• ,N - 1, 

.0 otherwise. 

The above currents can be either electric or magnetic the latter being 
regarded as equivalent to ideal electric voltage sources. 

Note that the spherical coordinate systems in Fig. 1 and 2 differ 
from those commonly used in phased array analysis. The poles are 
located at endfire instead of broadside and the ranges of colatitude and 
azimuth are such that the upper hemisphere is spanned by ^ 8 ^ r, 
^ <p < 7r. This convention is chosen for reasons of mathematical 
convenience. The results derived in Section III are valid for linear 
as well as planar arrays. The inclusion of both cases in a single 
treatment is facilitated by a generalized notation for the current 
density excitation function. The steering phases m\p and m\f/ x + 
nft are replaced by an equivalent "steering coefficient" a mn (<p vq ) in 
the plane of scan oriented at azimuth angle <p va . The steering coefficient 




Fig. 1 — Linear phased array. 
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Fig. 2 — Planar phased array. 

is derived by its relationship to the direction of a beam's maximum, 
which is determined for linear arrays by the equation 



p = 0, ±1, ±2, 



± 00 



(3) 



yp + 2pr — ka cos 6 pi) 
and for planar arrays by 

&, + 2p?r = A;a cos d VQ cos ^ p = 0, ±1, ±2, • • • ± oo (4a) 

ft, + 2qw - kb cos d PQ sin ^ PB g = 0, ±1, ±2, • • • ± », (4b) 

where A; is the wave number in the medium, and 8 pq is as shown in Fig. 
1 and 2. The steering coefficient is then defined by 

<r m n(<Pva) = k(ma cos <p pa + nb sin <p PQ ) 

p, q = 0, ±1, ±2, • •• ± oo. 
Equations (1) and (2) can now be written as 

J (a; — ma, y — nb, z) exp (j<r m „ cos d VQ ), 
ma ^ x ^ (m -\- \)a, 
nb g y ^ (n + 1)6, 
m = 0, 1, ••• , M - 1, 
n = 0,1, ■■■ ,N - 1, 
otherwise 



(5) 



T(as, ?/, 2, O = . 



(6) 



at <£ p( , = const. 
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Under the above generalization the excitation function for the 
linear array becomes a special case, q = 0, N = 1, <p VQ = 0, and the 
period in the {/-direction extends from — <*> to + °o ; or alternatively 
p = 0, M = 1, <p VQ = t/2 and the period in the ^-direction extending 
from -co to +oo. Since the phase constant exp \j(r m „(<p pa ) cos 6 pq \ 
is independent of (p, q), any p<1 may be chosen as the independent 
variable of scan. The subscript pq will henceforth be omitted whenever 
the mathematical expressions are independent of (p, q). 

The time dependence e' ut is assumed throughout the analysis. In 
a steerable array the phase taper is time dependent. However, it is 
understood that the rate of change of the phase taper is very small 
in comparison to the angular frequency, i.e., d$/dt « w, since only 
under that condition do the classical concepts of directivity and radia- 
tion impedance remain meaningful. If rf/(t) is a step function it is as- 
sumed that the time interval is long enough to allow all transients 
to reach negligible values before a new step is initiated. 

The formal solution of the array problem is obtained from Max- 
well's Equations via a vector potential A(x, y, z, 6) which is a solution 
of the inhomogeneous reduced wave equation 

V 2 A + k 2 A = -nj(x, y, z, d), (7) 

where n is the permeability of the medium. The magnetic field is given by 

H = - V X A, (8a) 

and the electric field (under Lorentz gauge) by 

E = -jJa + p VV-AJ- (8b) 

The solution to (7) over infinite space V can be written in closed 
form in terms of a dyadic Green's function 7 

A(.t, y,z,d) = n f S(x, y,z\S,y,l;)-J(t,V, t, 9) d£ d v d?, (9) 

where §(:r, y, z \ £, y, ?) is a solution of 

dx dy dz 

I being the unit dyadic a r a z + a.a,, + a 2 a, . The boundary conditions 
which § has to satisfy are derivable via the Vector Green's Theorem* 

* P. M. Morse and H. Feshbach, Ref. 7, p. 1767. 
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by imposition of the requirement that the tangential component of 
the electric field has to vanish on all conductors. This guarantees that 
all induced currents are accurately determined. 

It can be shown 8 that the average complex power delivered by the 
ninth, element in the array is 



P m „ = -= 



\j E-J* n dv, (11) 



where 

Jmn(x,y,z, 6) = J(x,y,z, d), 

ma£ x g (m + l)o, nb ^ y ^ (n + 1)6 (12) 

the asterisk (*) denotes complex conjugate, and V mn is a simply con- 
nected volume occupied by J m „ . If S mn is a surface obtained by taking 
a cross section through V mn , the total current, I mn , flowing through 
the cross section S mn is 



-// J-ds. (13) 



The element radiation impedance, Z mn , is denned in terms of the com- 
plex power by 

P = * | I mn | 2 Z mn . (14) 

By (10) and (13) via (8b) and (9), the element radiation impedance 
can be defined directly in terms of the array geometry and the excita- 
tion: 

zUe) = TT—r- f f J5-<*« •''• 2 - *>'&(*' v> » I *. i. A 

| 'mn | J Vmn J V 

■1,1, f, e)drdv, (15a) 
where dr = d£<hjd(, dv = dxdydz, and 

G(x, y, z | f, ,, t) = m\I + p w)-9(.r, s,, z | $ , „, f). (15b) 

Operator V operates on (.r, 7/, 2). The quantity |/ m „| 2 is introduced for 
the purpose of normalization, and may depend on the choice of the 
cross section Smn. 

The definition of the impedance includes both linear and planar ar- 
ray elements. It is consistent with the commonly known definition of 
impedance" if the latter is viewed as a relation between the average 
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complex power delivered by the generator and the rms current flowing 
into the load. The definition given by (15) is necessary in view of the 
fact that in a system excited by distributed currents, a terminal volt- 
age in the time domain is not always uniquely defined. In a system 
excited by magnetic currents, (15) defines the element admittance if 
the permeability /j, is replaced by the permittivity e and the electric 
currents by their magnetic counterparts. 

In the following theoretical discussion, it is assumed that the phased 
arrays are excited by a uniform amplitude and a linear phase taper. 

III. FINITE ARRAYS 

Theorem 1 : The element radiation impedance in a finite, steerable, linear 
or planar phased array of scan-invariant current sources, radiating into 
a linear, lossless, passive and time-invariant system, is an entire function 1 
of the scan angle 6 in any given plane of scan, with an essential singidarity* 

at 6 — > co . 

Proof: By (15a) 

I Jr. \ 2 Z rs = f f Jr*.(*. V, 2, d)-Q(x,y,Z | *, V , f) ■ J(£, t, , f - *) *T dv . 

(16) 

On expanding (16) in a double sum of integrals over all cells { (m, n) }, 
m = 0, 1, • • •, M — 1; n = 0, 1, • • •, N — 1, and using the relationships 
of (6) followed by a change of variable in each term of the sum, one 
obtains 

Zrs(d) = Z IU exp [j(<r mn - O cos 6], (17a) 



• n=0 



where 

2m„r» = 



f f Jt(x,y,z) 



■G(x + ra, y + sb,z\S+ ma, y + nb, f) 'J ft, n, f) dr dv. (17b) 
In any given plane of scan <p is constant, so that 

o-„„, — <x rs = k[(m — r)a cos <p + (n — s)6sin<^] = <r m - r , n - t (18) 
is independent of 6. Both cos $ and the exponential function are entire 
functions. 1 " Consequently, the exponential function appearing in (17a) 

* R. V. Churchill, Ref. 10, Sec. 68, p. 157; Sec. 112, p. 270. 
t R. V. Churchill, Ref. 10, Sec. 21, p. 47; Sec. 23, p. 50. 
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is an entire function of an entire function, which is likewise entire 11 
(entire functions are also called integral functions). Z ra (8) is a finite 
sum of entire functions and is also entire. 

The nature of the essential singularity at -» oo is obtained by first 
expanding cos fl in the complex 0-plane 

cos (0 r + jd.) = cos r cosh 0, - j sin r sinh 0, . (19) 

Then, if | t | — > <x in such a way that {c mn — a r .)di > 0, the ranth 
term behaves as exp {[ o mn — o r , | sin r exp [| 0, |]} Q.E.D. Note 
that even when J (x, y, z, d) is scan dependent, Z r ,(d) is analytic pro- 
vided J (.r, y, z, 0) is analytic. However, other isolated singularities 
may exist. 

Corollary la: Re \Z r ,\ and — Re{Z r ,\ are entire junctions of each 

da 

with an essential singularity at d — > x . Proof appears in Appendix A. 

Theorem 2: The power radiated by an element in a finite, steerable, linear 
or planar phased array of scan-invariant current sources, radiating into 
a lossless, linear, passive and time-invariant system cannot be kept con- 
stant over a continuous scanning range with lossless, linear, passive and 
time-invariant network elements and scalterers only. 

Proof: Let G(x, y, z j & rj, f) be the dyadic Green's function of the 
entire system including all equalizing elements. The radiation impedance 
of the wmth element of the array is given by (15a) for a lossless, linear, 
passive, time-invariant system. If the array is radiating constant power 
over a continuous scanning range, the real part of the radiation im- 
pedance, R r ,{6) = Re{Z r . }, must remain constant in that range and 

— [R ri {0)\ = o, 0, ^ e r ^ e 2 , e, ■ = o (20) 

where 6 = 6 r + J0, . By Corollary la, -^ [R„(d)] is analytic in the 

(70 

closed 0-plane and has an essential singularity at — > oo . However, 
if the derivative vanishes along the line 0, ^ r 5S 2 it must vanish 
everywhere in the 0-planc*. Hence, it cannot have an essential singu- 
larity at infinity. The contradiction implies that R r ,(6) cannot be 
constant over a continuous scanning range. Q.E.D. 

Equations (3) and (4) specify the directions of the beams' maxima, 
however, not all of them correspond to real directions in space. Whereas 
<fpq is real for all {p, c/) , 6 m can be either real or imaginary, as may be 

* P. M. Mor.se and H. Feshbach, Rpf. 7, Vol. I, p. 390. 
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seen from the solution of (4) : 

*.-^ft£8i- °**- < - (21a) 

If pa is real it is said that the beam is in real space. By way of mathe- 
matical generalization it is said that all those beams having an imaginary 
$„ are in "imaginary space". If P , = 0, or d p „ = ir, it is said that the 
beam is in a grazing position between real and imaginary space. It 
can easily be verified from (21) that for a given phasing (\p x , \p v ) every 
pair (p, q) corresponds to a unique direction (<p PQ , O in the complex 
domain ^ < f, ^ Re(fl| £ ir. These directions are the char- 
acteristic directions of the system. They are directly related, through (4), 
to the eigenvalues of 

Q + S + T2f ( x > *> - ° (22) 

dx dy 

with the following periodic boundary conditions 

F(x, y) = F(x + a,y) exp (-;>■), ( 23a ) 

fj (.t, y) = ^ (x + a, y) exp (-;>.) , (23b) 

F(x, y) = F(x, y + b) exp (~M, (23c) 

^ (.r, y) - ^ (*, y + 6) exp (-;>,). (23d) 

The eigenfunctions, which form a complete orthogonal set in the 
interval ^ x ^ a, ^ y ^ & are 

IWr, 2/) = exp [;(*, + 2?*r) |] exp [;(*„ + 2r 7 7r) |J , 

p, 5 = 0, ±1, ±2, ••• ± oo. (24) 
By (4) they can also be written as 

F va (x, y) = exp [jk cos fl^x cos <p vq + y sin $>„)}. (25) 

The eigenvalues {r pg } are 

r DB = k cos e va p, Q = 0, ±1, ±2, • • • ± co . (26) 
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The results thus derived lead to several interesting conclusions which 
are summarized in the following lemmas. 

Lemma 1: Every steerable linear or -planar phased array with a linear 
phase taper has only a finite number of beams in real space. Proof appears 
in Appendix B. 

For every pair of phasing (\f/ z , \p v ) there exists an infinite set of 
characteristic directions [d PQ , <p PQ }. As the array is scanned by varying 
the values of (\p x , \f/ y ) in the intervals — ir ^ \f/ a ^ r, — ir ^ \J/ V ^ ir 
some characteristic directions will go through a grazing position going 
from imaginary to real space or vice versa. We shall call such char- 
acteristic directions "transitive characteristic directions".* Since the 
condition for a grazing position is | cos 6 PQ \ = 1, it follows from Lemma 1 
that the number of transitive characteristic directions is finite. 

Lemma 2: The radiation impedance of an element in a linear or planar 
phased array can be expanded by an infinite series over all characteristic 
directions of the system. Proof appears in Appendix C. 

IV. INFINITE ARRAYS 

In analyzing large arrays it has been found useful to approximate 
the behavior of the center elements by the behavior of identical ele- 
ments in an infinite array of the same geometry. 12 This approximation 
is motivated by the fact that the performance of the center elements 
is strongly affected through mutual coupling by contiguous elements, 
but very weakly by elements far away. 13 

The formulation of the infinite array problem may be obtained from 
the results derived for finite-size arrays by letting the number of 
elements M and N approach infinity. The infinite array problem can 
also be treated as a periodic structure by application of Floquet's 
theorem. In the following, the latter approach is adopted, but first it 
is demonstrated that both methods are consistent. 

The electric field of an infinite array as given by (8b) must satisfy 
the same periodicity conditions as the source function, i.e., 

E(.r + ma, y + nb, z) = E(.t, y, z) exp [j<r m „(<p) cos 6). (27) 



* Note the distinction made between "grazing position" and "transitive charac- 
teristic direction". A beam associated with a transitive characteristic direction 
may attain a grazing position for a particular phasing, but may also point in other 
directions. 
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On the other hand, the electric field 

E(a,-, y, z, 6) - - f G(x, y, z\£,v, f)-J(*. V, f, 6) dr (28) 

can be expanded in an infinite sura of integrals using the relationships 
of (6) : 

00 00 

E(x, y, z, 0) = - E E exp (j<r m „ cos 6) 

■ f G(x, y,z\Z+ ma, v + n&, f) • J ft, *, f) *, (29) 

where F o is the volume occupied by J . Define a new Green's function 
(x, y,z | £, v, f) 

= E E exp 0V mn cos 0)<5(.r, y, « | £ + «m, 77 + nb, f) (30) 

m = — 00 n — — 00 

and notice that 

&(*. !/,2|? + Ma, v + #&, f) = exp (—faun cos 0)&(s, ft z | £, 17, f) 

(31) 

since by (5) 

O m + M,n + N = <*mn + ^JtfAT • ( 32 ) 

From (27) and (31) it follows that <5 (a;, ft « | £, 17, f) can be expanded 

by the eigenfunctions (25) as 

&- E E ^ flJJ* id«<l. ^. ( 33a ) 

JJ--00 Q--W 

where 

fc, - 5 £ £ 8oF»fo v)F*(x, y) dx dy dt- d v . (33b) 

Substituting (30) via (33a) into (15a) for the center element, m = 
n = 0, one obtains 

Zoo = E E *„ , ( 34a ) 

p=— 00 a-— 00 

where 

z pa = 77 — [5 / / J?(«, 2/. 2) 

I JOO I •'T'oo •'Voo 

• Sp.fr I f)/%(€, l)*U*. y) ' Jo(l, r?, f) dr cfo. (34b) 
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Equation (34) is an alternate representation to (86) for the radiation 
impedance of the infinite array element and it demonstrates that 
Lemma 2 is valid for infinite arrays as well. 

By substituting the new representation for G , (30), (33), into (29) 
and noting that the electric field satisfies the homogeneous reduced 
wave equation in the source-free region, one obtains for the unbounded 
space 

oo co 

E(s, y, z) = Z Z SpAafr, y) exp (-y pq \z\) \z\> d mnx , (35a) 

j»— — oo a-— oo 

where 

7 P , = Vr 2 pq - P = jk sin e pv (35b) 

S M <f 1, " ,,l = - f g„(z | fl.Joft, *, iWfc, v) rfr, (35c) 

and d max is the projection on the z-axis of the largest distance between 
two points on the surface enclosing F o- It can be seen that the electric 
field in the source-free region, above the central area of a large array 
may be approximated by a finite number of homogeneous plane waves 
propagating in the real characteristic directions, and an infinite num- 
ber of nonhomogeneous plane waves, exhibiting exponential decay in 
the direction perpendicular to the plane of the array. The latter are 
interpreted as waves propagating in the imaginary characteristic di- 
rections. 

In an infinite array all elements are embedded in an identical en- 
vironment, and therefore the power radiated by each element is the 
same. There is no net power flow into a unit cell through the "side 
walls". Consequently, the quantity Re{|Joo| 2 2p 9 } of (34b) is equal to 
the power propagated by the plane wave (p, q) within a unit cell in 
the direction perpendicular to the plane of the array. By Lemma 1 
there is only a finite number of plane waves with transitive charac- 
teristic directions (see footnote p. 1571). Let them be distinguished 
from all other plane waves by assignment of the subscript (p, q) = 
(t,v). 

E r , = B„F t ,(x, y) exp (-jk | z \ sin B„), \z\> <4.x (36) 

H T „ = 3C Ty F„(x, y) exp (-jk\z\ sin 6„), | z \ > d mKX , (37) 

where F„(x, y) is given by (26), and £ r „ by (35c). If 

S>„ = jk[co3 6 T , cos(p Tt a. x + cos 6 Tr sin <p Tt a„ — sin T „aJ (38) 
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then 

3C r , = X ©„ X 8„ . (39) 

OljU 

The power radiated by a (t, v) plane wave per unit cell into the upper 
hemisphere is 

P T¥ = hie f f (E„ X H*).a, <fe <Zy. (40) 

2 Jo Jo 

Substitution of (36) through (39) into (40) gives 



n ab . * 

P rp = — sin 9„ 

Vo 



,9 , . 12 i sin T „ , _ 



(41) 
where ^ = W«)*- From (41) a radiation resistance per wave is defined as 



R„ = 7^72 • (42) 

Moo I 

Since the entire system is passive and lossless, then by conservation 
of energy, the power P TV must originate from the element itself. Hence, 

R„ = Re {z„), 0„ real, (43) 

where z r „ is given by (34b). 

From (41) it follows that when a wave (r, v) is in real space R„ 
is real, and when it is in imaginary space R r , is imaginary (in which 
case Re{z T „) = 0). Hence, of all the elements comprising the source's 
load, R Tr appears either resistive or reactive, depending upon the 
scan angle. Such properties of a load, which are unknown in lumped 
network theory, are a consequence of the losslessness postulate. When 
propagation is possible power is carried away from the source. When 
propagation is inhibited there is no net loss of power and the load 
must be reactive. By Lemma 1 only Re{z r ,\ has those properties. 
All other z pg , (/>, q) ^ (r, v) and Im{z„} always retain their dis- 
sipative or reactive characteristics. Further, there is only a finite 
number of terms having Re{z pa } > 0. In practical phased arrays the 
spacing between the elements and the scanning range are such that 
only one such term exists at a time. 

The following two definitions summarize the properties described 
above: 

Definition 1 : An O-type network function is a scan-dependent immit- 
tance (impedance or admittance) which is seen by the source as resis- 
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tive when the beam is in real space and as reactive when the beam is 
in imaginary space, and it behaves like an open circuit for impedance 
and like a short circuit for admittance in the grazing position. 

Definition 2: An E-type network function is a scan-dependent immit- 
tance (impedance or admittance) which remains either resistive or re- 
active when the beam passes through the grazing position. 

The motivation behind the nomenclature introduced by the two defi- 
nitions will become clear later, in Theorems 4 and 5. The O-type and 
E-type immittances are of distinct mathematical form. To arrive at it 
consider first the following transformation:* 

s = sin 6 mn (44) 

X = cos <p mn , (45) 

where (m, n) is one particular transitive characteristic direction out 
of all (t, v). Given s and x a ^ other characteristic directions are 
uniquely determined. By (4) 

.£ x = kaxVl - s 2 - 2mw (46) 



yp v = kbVl - x Vl - s" - 2mr, (47) 

where (1 - x 2 )' ^ for all possible x and (1 - s 2 )* > Oif ^ 6 mn < tt/2, 
and (1 - s 2 ) h < if tt/2 < 6 mn ^ t. Then by substitution of (47) 
into (22) all other characteristic directions are found: 

tan fc _ ^ a-fa-^ + ^-n)™ 

kabxO- — sy + 2(p — m)irb 



cos $ pq = f„(s), (48b) 



where 



ka cos <p PQ K ' 

This suggests that when characteristic direction (m, n) is scanned in a 
plane x = const, each of the components z PQ of the total input impedance 
as given by (34) can be expressed as a function of the same variable s. 
The conformal mapping between the m „-plane and the s-plane is shown 
in Fig. 3. In view of the branch cut — 1 ^ a ^ 1 it will be understood 
that s = a denotes s = a — ;0 if ^ 9 mn ^ tt/2 and s = a 4- JO if 
ir/2 ^ 0,„„ ^ ir. Let s = s T , be the value at which characteristic direc- 
t-Recall that 0„i >i and (p mn arc not in the conventional spherical coordinate sys- 
tem (see p. 1564). 
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lm{0 mn } 

c 




Fig. 3 — Conformal mapping s = sin mn . 

tion (t, v) is in grazing position. At this value 

f„(8„) = 1. (49) 

Of all values {s T „\ there is at least one which satisfies (49) for s r , = 0. 
From (48) it is obvious that f mn (0) = 1, and there may be other transi- 
tive characteristic directions (r, v) 9± (in, n) which attain their grazing 
positions at s T , = 0. 

Theorem 3: In an obstacle-free space, the impedance junction z PQ (s), 
associated with the characteristic direction (p, q), is an analytic function 
of the complex variable s = a + j0, with branch points at s = s T „ and an 
essential singularity at \ 8 | ->«>.// (p, q) = (m, n) then z nn (s) may 
have a simple pole at s = s„,„ = 0. 

Proof: The general definition of z PQ is given by (34b) in which the 
pa dependence is contained in the Green's function component 
-g pq {z | f)F p *(£, v)F PQ (x, y). The Green's function is derived from (10) 
via (15b). Green's function §(x, y, z \ £, v, f) satisfies the same periodic 
boundary conditions as G Q (x, y, z | £, v, f) and can be expanded in a 
series similar to (33a) : 



9= Z Z C pa (z\^F va (x,y)F:^, v ). 



(50) 



j> =— oo a - — oo 
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By substitution of (50) into (10) and use of the orthogonality property 
of F VQ {x, y) one obtains a differential equation for C pv (z | f) 

d? -T-^ltf - - J — -ab~^ 

(51) 

y PQ = jk sin P , 

with the additional requirement that as | z | — > °°, (7 P , behaves as an 
outgoing or evanescent wave. The solution of (51) for free space is 

" (Z I f) " ' ijabk sin , m CXP ( ~ ik I Z - f I Sin "" I < 52 > 

<7 P «(2 | f) is obtained from C p „(z | f) through an operator (R„, : 

5„,(2 | f) = pH^-CJ? | f), (53a) 



rhere 



(Rp„ - / + p »*„©,,„ (53b) 



S) p0 being given by (38). Substitution of (52) into (53a) followed by 
substitution into (34b) gives 



S I lT\ 2 Sin eZ Jy„ ]y„ JUX ' V ' 



1 -2a hi ■ 

'<Rm-Jo(£i V, f) exp {;'& cos 6 pq [(x - £) cos «p P9 

+ (2/ - i?) sin <p„] - jk | 2 - f | sin 0,.} dr dw. (54) 

The integrand is an entire function of d pq with an essential singularity 
at | Im \0 PQ \ | — > oo. Hence,* if J (.r, ;/, 2) is piecewise continuous, 
the integral is also an entire function with the same essential singularity. 
By (48) 

sin 0„ - [1 - &(*)]*. (55) 

By Lemma 1, 

&(•) 3* 1 if (p,<?) * (r.r). (5(5) 

From Fig. 3 it is readily seen that | 8 | < 00 when | d UQ \ < 00 which 
implies, via (48), (55) that | cos P « | < °° and | sin P(Z | < co as long 
as I s I < co . Thus, the singularities introduced by the transformation 
(44), (45) are the branch points at s = s„ . Also if (p, q) = (m, n)f 

* E. J. Copson, Ref. 11, Sec. 5.5, pp. 107-109. 

t Recall that (m, n) is the characteristic direction which defines the transforma- 
tion from (^«,^») into (s. x), (44)-(47). 
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z p9 may have a simple pole at s = 0. (Note, for example, that for hori- 
zontal polarization, J = a. x J , there is a simple zero in the plane of 
scan corresponding to <p mn = 0, at s = 0.) Q.E.D. 

The above proof can be applied separately to the real and imaginary 
parts of the right-hand side of (54). If z PQ = R ptt (s) + jX PQ (s), then R PQ (s) 
and X pv (s) are analytic functions of s, real on the real axis of s, with 
an essential singularity at | s | — > °o , branch points at s = s Tr , and 
possibly simple poles at s = 0. 

In systems other than obstacle-free space, the normalized complex 
power z pg (s) has different forms. Except for isolated values of s, the 
radiated power and the stored energy per unit cell are bounded and 
continuous functions of s over those portions of the real and imaginary 
axes of the s-plane which have physical meaning. Hence, it is reasonable 
to postulate that an analytic continuation of z pq as a function of scan 
can be made into a region of the complex s-plane which includes por- 
tions of both the real and imaginary axes. It may be of interest to 
note that the impedance function z pq (s) derived by L. Stark 14 for the 
planar dipole array over a ground plane is analytic. The regularity 
of z vq (s) depends directly on the regularity of g PQ (z | f ; s). The singu- 
larities of z pq in the s-plane are determined by the boundary conditions 
which g pq (z | f ; s) satisfies. 

Theorem 4: An E-type immittance junction V(s) is an even junction of s. 

Prooj: Let the complex variable s be defined with respect to the transi- 
tive characteristic direction (to, n). Once (to, n) is chosen, the proper 
branch of (1 - s 2 )* in (48) is uniquely determined. Let (k, I) denote 
all other transitive characteristic directions which reach their transitive 
position simultaneously with (?n, n). Formally, this implies 

/;,(0) = 1 (r, v) = (m,n), (k,l). (57) 

As a consequence of Definition 2 and Lemma 1, V(s) is recognizable as 

T v s) = [#*» (P> ?) * ( m ' n )' &< ®> (58) 



[X pq (s) all (p,g), 
where R pq (s) + jX pq (s) = z pq (s),z pq given by (54). Thus, (58) estab- 
lishes the connection between the defined E-type function and physical 
quantities corresponding to R PQ (s) and X pq {s). Consider Definition 2 
which states 

V(s) - V*(s) = s = a 0<a<l, (59a) 

V(s) - V*(s) =0 s = ;'/3. (59b) 
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Since V{s) is analytic and also real on the real axis of s, (59) may be 
rewritten as* 

V(s) - V(s*) =0 s = a < a < 1, (60a) 

Y(s) - Y(s*) = s = ft. (60b) 

On the real axis 

V(a) - Y(a) = 0. (61a) 

On the imaginary axis 

Vm - V(-j$) = 0. (61b) 

By analytic continuation 1 of (61b) from the imaginary axis to a point s 
in the complex plane one obtains 

V(s) - V(-s) = 0. (62) 

Hence, V(s) is an even function of s. Q.E.D. 

Theorem 5: An O-type immittance W(s) is an odd junction of s. The 
proof is similar to that of Theorem 4 and it appears in Appendix D. 

It has been shown in Theorem 2 that a finite phased array cannot 
be perfectly matched over a continuous scanning range. The proof is 
limited to finite arrays and cannot be directly extended to infinite ar- 
rays since the representation of the element impedance by (17a) does 
not guarantee convergence in the complex 0-plane if the limits of the 
summations are extended to infinity. In treating the infinite array, the 
element impedance is derived by symmetry considerations from which 
it is concluded that the net complex power radiated from each element 
is conserved entirely within the unit cell of that element. It has been 
shown that the two definitions are consistent. Although the problem of 
whether an infinite array can be perfectly matched is of academic in- 
terest only, it is worthwhile noting that as for finite arrays, the answer 
in this case is also negative. To show this the reader may recall that 
the impedance has been defined as normalized power and postulated to 
be an analytic function of the scan variable s = «+ jfi. The normaliza- 
tion constant is |/ o| 2 given by (13). If the complex power as a func- 
tion of scan is represented by 

Pis) = \I 00 \ 2 [R(s) + jX( S )], (63) 



* P. M. Morse and H. Feshbaeh, Rcf. 7, Vol. I, p. 393. 
t Morse and Feshbaeh, Op. Cit., p. 389. 
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then by Lemma 1, the term R(s) is a finite sum of analytic functions 
of the complex variable s. Consequently, R(s) is an analytic function 
of s. In general, it may be represented as 

R(s) = E(s) +0(8), (64) 

where E(s) is an even function of s and 0(s) is an odd function of s. 
Under conditions of perfect match over a continuous range, constant 
power, P r , is radiated over that range. Since R(s) is analytic it implies 
R( s ) = p r | 7 00 |~ 2 everywhere in the s-plane. Since a constant is even, 
0(s) = 0. Further, E(s) must have a branch cut on the real axis of the 
s-plane in the interval [—1,1]. But the branch cut does not exist if 
E(s) = P r | Zoo | _i . The contradiction implies that P(s) in (63) cannot 
equal a constant over a continuous range of s. 

Theorem 6: The resistance and reactance junctions oj an element, or their 
derivatives, in an infinite linear or planar -phased array of current sources 
are discontinuous when a grating lobe is in a grazing position. 

Proof: In an infinite array the grating lobes are plane waves propagating 
in the characteristic directions. By Theorems 4 and 5 the element 
impedance Z(s) in an obstacle-free space can be written as 

Z(s) = P(s) + 2®- (65) 

For real values of s, P{s) is an even complex function of s bounded at 
s = 0, and Q(s) is an even real function of s nonzero at 8 = 0. On the 
real axis of s 

Z(a) = P(a) + ®&- (66a) 

On the imaginary axis of s 

zm=pm-j g f l - (66b) 

A grating lobe is in its transitive position at s = 0. The pole discon- 
tinuities are established by showing that 

Re {lim Z(a) - lim Z(j/3)} = lim ^ = °o (67a) 

a-,0 0-0 a-0 OC 

Im (limZ(a) - lim Z(j$)\ = lim ^^ = °o . (67b) 

O-.0 0-«O p-0 P 
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The pole discontinuity has to be interpreted as an invalid mathe- 
matical solution at the transitive position. It is a result of the idealiza- 
tion introduced by the concept of an "infinite array." If R mn (s) has a 
simple zero at s = 0, as is the case when a horizontally polarized array 
is placed above a ground plane, then the active impedance in the 
neighborhood of 8 = can be written as 

Z(s) =R(s) +jX(s), (68a) 

where B{s) and X(s) are real functions of s (real for s real). 

R( s ) = f) a #* (68b) 

X(s) = E b 2i s 2i . (68c) 

• ■=0 

When the beam whose transitive characteristic direction is in real 
space, s = a 

R a 4 /f(a ) = £ a . a i ( 69a ) 

i=0 

and when it is in imaginary space, s = j/3 

R fi ± Re \R(j(3)} = E (-l)W'- (69b) 

i=0 

The discontinuity in the derivative of the resistance is 

hm —j hm -rf = a x . (70) 

a _ da p^o dp 

Similarly, the reactance 

X a ± X(a) = E 6 2 ,a 2 '' (71a) 

> =n 

X, = Ira [Z(j0)} = E (-D'lW + fl 2 . + ^ 2,+1 ] (71b) 

i -0 

. . dX n . . dXg /» \ 

hm— lim-rr- = — a, . \Jz) 

a _.o a a ff^n up 

The proof can be generalized for any order algebraic singularity or 
zero at s = 0. For example, if there is a zero of multiplicity N the dis- 
continuity will be in the Nth derivatives of the resistance and react- 
ance. A noninteger order zero yields a discontinuity after a sufficient 
number of differentiations. Q.RD. 
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V. SUMMARY AND CONCLUSIONS 

A new mathematical approach to phased arrays has been adopted 
to investigate and discover various properties of the radiation imped- 
ance of an array element as a function of scan angle. The underlying 
idea of the method is the treatment of the impedance as an analytic 
function of a complex scan variable, which enables one to prove that 
an array element subject to the model chosen cannot be perfectly 
matched over a continuous scanning range by using lossless, linear, 
passive and time-invariant elements. 

The first half of the theory is devoted to finite arrays. It is shown 
that the directions (in space) of the beams' maxima are eigenvalues 
of a Laplacian differential operator with periodic boundary conditions, 
which are related to the phase taper of the array. It is proven that 
there exists only a finite number of real eigenvalues. The known con- 
cept of imaginary space is then adopted to accommodate the imagi- 
nary eigenvalues. Furthermore, it is demonstrated that all beams except 
a finite number are completely confined either to real space or to 
imaginary space, and that only a finite number of beams may attain 
a grazing position. The unique properties of the latter beams have 
been found to play an important role in the investigation of infinite 
arrays, to which the second half of the theory is devoted. 

The interest in infinite arrays, apart from its academic aspect, stems 
from the good approximation it provides for the behavior of the cen- 
ter portion of a large finite array. It has been found that the infinite 
array element impedance as a function of scan is restricted to a spe- 
cific mathematical form. It is the authors' hope that recognition of the 
limitations imposed by that form may provide useful guidelines in 
achieving optimal match of an array to space. 
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APPENDIX A 

Proof of Corollary la 

Corollary la: Re{Z T ,\ and — Re{Z r ,\ are entire functions of 6 each 

ov 

with an essential singularity at d —> «> . 
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Proof: Denoting 



one obtains from (17a) 

Re{Z ri (0)} ±R r .(e) 
jvr-i n-i 

= EE (p. 



where 



and 



B, 



2mnr. = P*nr. + jXmnr, 

„, cos [£ mnr .(ff)] - x„„,. sin [B n 

mr.(^) = (ff«» — O COS 5 



.(^)] 



(73) 

(74) 
(75) 



Off m-0 n=U 

+ x cos[i? mnrB (^)]!. (76) 

Since cos 9 is an entire function of 6, cos[B„, nrs {$) ] and sin [Z? „„„..„ (0) ] 
are entire functions of an entire function, and are therefore entire. The 
existence of the essential singularity can be demonstrated in a similar 
fashion to that in Theorem 1. Q.E.D. 
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Proof of Lemma 1 

Lemma 1 : Every steerable linear or planar phased array with a linear 

phase taper has only a finite number of beams in real space. 

Proof: A beam (p, q) is in real space if | cos M | ^ 1. Dividing (4a) by 
ka and (4b) by kb, squaring and adding, one obtains 



ka I " r \ kb 



^ 1 



or 



Necessary conditions for the above inequality to be satisfied are 

a 



2tt 


+ 


V 


2tt 


+ 


S 



(77) 
(78) 

(79) 
(80) 
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Since 

2 - 2tt = 2 an(J 2 = 2tt = 2 ' 

bl^f + | (81) 

l*l*X + 5" (82) 

Hence, both 7) and <7 arc bounded. Q.E.D. 

appendix c 

Proo/ of Lemma 2 

Lemma 2: The radiation impedance of an element in a linear or planar 
phased array can be expanded by an infinite series over all character- 
istic directions of the system. 

Proof: The current density excitation function of a finite-size array 
given by (1), (2) satisfies the periodic boundary conditions (23) in 
the finite domain occupied by the array. Let this domain be denoted 
by D. The current density can, therefore, be uniquely expanded in D 
in terms of the eigenfunctions (25) : 

J(x,y,z) = U(x,y,D) £ fj j„(*)P M (a:, y), (83) 

where 

)M = -T [ " f W*. V. *W«(*. tf * * (84) 

£10 Jo ^o 

and £/(z, ?/, D) is a two-dimensional unit step function 

t/(«. Vl D) = j 1 (*■»> W D ' (85) 

[0 otherwise. 

Substitution of (31a) into (15a) yields 

Z mn = t, Z «.«. . (86a) 

where 

•J w (f)«l)^.U.O)^*. (86b) 
Q.E.D. 
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APPENDIX D 



Proof of Theorem 5 

Theorem 5: An O-type immittancc W(s) is an odd function of s. 

Proof: Let the complex variable s be defined with respect to the tran- 
sitive characteristic direction {m, n) . Let (k, I) be all other transitive 
characteristic directions which reach their transitive position simul- 
taneously with (m, n). Then as a consequence of Definition 1 and 
Lemma 1 

W(s) = R P M (p, q) = (m, n), (k, I), (87) 

where R pg (s) = Re{z„ q }, z pg given by (54). Thus, (87) establishes the 
connection between the defined O-type function and a physical quan- 
tity corresponding to R PQ {s). From Definition 1 

W(s) - W*(s) = s = a, 0<«<1 (88) 

W(s) + W*(s) = s = i/3. (89) 

Since W(s) is real on the real axis of s, (88), (89) may be rewritten as 

W(s) - W(s*) = s = a (90) 

W(s) + W(s*) = s = ;/3. (91) 

On the real axis 

W(a) - W(a) = 0. (92) 

On the imaginary axis 

W(jp) + !F(-i/3) = 0. (93) 

By analytic continuation of (93) from the imaginary axis to a point 
s in the complex plane one obtains 

W(s) + TF(-s) = 0. (94) 

Hence, W (s) is an odd function of s. Q.E.D. 
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